
3 1 


76 01404 0928 


UA5A-TM-?ms 


NASA TECHNICAL MEMORANDUM 


NASA TM- 88445 


NAS A-TM-8 8445 

ITERATIVE ADAPTION OF THE BIDIMENS I ONAL WALL OF THE 
FRENCH T2 WIND TUNNEL AROUND A C5 AX I SYMMETRICAL 
MODEL: INFINITE VARIATION OF THE MACH NUMBER AT 
ZERO INCIDENCE AND A TEST AT INCREASED INCIDENCE' 


J.P. Archambaud, J . B , Dor, M.J. Payry, 
L . Lamarche 


Translation of "Adaptation iterative des parois 
bidimensionnelles de la soufflerie T2 autour 
du modele axisymetrique C5 : variation du nombre 
de Mach infini a incidence nulle et un essai de 
mise en incidence". Rapport Technique OA 35/ 
3075 (DERAT No. 13/5015) DNO, O.N.E.R.A. , Centre 
d' Etudes • et de Recherches de Toulouse, France, 
March 1986, pp 1 - 16 (plus figures) 



JUL 2 1986 

LANGLEY RESEARCH CENTER 
LIBRARY, NASA 
HAMPTON, VIRGINIA! 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON D.C. 20546 JUNE 1986 


/ 



NF00983 


I 


stawoa#© nvLg PA(jg 


5 W©*- l ©-©irowfiiwsf alert Mo* 

i/fsstr iy\ -■SW! 1 __ _ _ 

3, <2 o c 1 pi 9vif * o Co*cto$ N®. | 

«. Ytile S«bmlo ITERATIVE ADAPTION OP THE BI- . 

DIMENSIONAL WALL OF THE FRENCH T2 WIND TUNNEL 
AROUND A C5 AXISYMMETRICAL MODEL, 

3. ftoperf D q Jo | 

JUNE 1986 

G. P •rfomilng Orgonll«H»* Cod© 

7 , Au*fto'(o) 

J.P. Archambaud, J.B. Dor, M.J. Payry, 
L. Lamarche 

0. pBff*'n»mg Orgoni notion M©, 

10. Wo,b Unit No. 

H, ContrgCP Of Grenf N*9» 

HASw 40 0 4 

y. Orgorii €/>d A 

SCITRAN * 

Box 5456 

Santa Ba.rha.ra, t'A 9 11 OR . . _ 

1$. T yp© of U v p © ? » «i)(j Pc»i©<2 Cov®r©<fl 

TraaoXation 

12. Sy i»f* son ng Aoonty AdJfOOO 

NLiitionai Aerocuiutieu cad Spaoa Adtiiniatsratioa 
I Wa8nlii3Con„ D.C. 20546 


50. t^onaevir.Q Agcvsjy C«nJ® 




!$, Net^a 


Translation of "Adaptation iterative des parois bidimension- 
nelles de la soufflerie T2 autour du modele axisymetrique 
C5" variation du nombre de Mach infini a incidence nulle et 
un essai de mise en incidence". Rapport Technique OA 35/ 

3075 (DERAT No. 13/5015 DNO, O.N.E.R.A., Centre d'Etudes et 
de Recherches de Touloqse, France, March 1986, pp 1 - 16 < 

]&, Ab«»r©«t' 

The top and bottom two-dimensional walls of the T2 wind 
tunnel are adapted through an iterative process. The 
adaptation calculation takes into account the flow three- 
dimensionality. This method makes it possible to 
start with any shape of walls. The tests were carried 
out with a C5 axisymmetric model at ambient tempera- 
ture. Comparisons are made with the results of a 
true three-dimensional adaptation. 
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Author's summary: During these tests the top and bottom two-dimensional 
walls of the T2 wind tunnel are adapted through iterative process. The 
adaptation calculation takes into account the flow's _thre^~dimensionali ty ; 
upoii each repetition the interference of the four walls at the model is 
* Numbers in. margin refer to foreign pagination. 
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reduced. This method therefore enables to start with any shape of walls. 
Another advantage is that the walls are adapted around a bearing model. 

The tests were carried out with a C5 axisymmetric model (braking 
coefficient 1.8%, length 400 mm) at ambient temperature and an average 
impact pressure of 1.7 bars (Reynolds number based on the length of the 
body: 10? at Mo=0.84). They include: 

- A Mach number sweep from 0.6 to 0.95 with zero incidence 
•- An adaptation at incidence (X= 8° at Mo=0,84 

- A non-converged test with the model greatly off center with 
zero incidence. 

For certain Mo values with zero incidence comparisons are made with the 
results of a true three-dimensional adaptation (TU Berlin) and tests 
conducted in a large wind tunnel (NASA Ames). 
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1 - INTRODUCTION 

■i An initial test series (September 1984, ref. 1) had shown the validity 
of two-dimensional adaptation around three-dimensional models. The 
adaptation method was based on a speed correction calculated in one step in 
three dimensions. 

The present series is a continuation of the preceding one and extends 
the range of validity of the process used thanks to the following 
properties: 

- The obstacle considered is the C5 axisymmetric body (ref. 2), already 
tested in 1984. But this time the model is larger (C - 400 mm; braking 
coefficient 1 . 8 %). 

- The calculation method developed by L. Lamarche (ref. 3) differs from 
the first version used: the method is a more precise iterative one which 
enables the interferences of the walls surrounding the model in incidence 
or off center to be corrected. 

2 - WIND TUNNEL T2 AND MODEL C5 

Wind tunnel T2 is an induction-driven closed system cryogenic facility 
with winds of 1 to 2 minutes long (fig. l) - Refs. 4 and 5. The present 
testing series was conducted only at ambient temperature. 

The test channel (H = 370 mm, B = 390 mrn at the inlet) comprises two 
flexible walls, one on the top and one on the bottom; each wall's shape can 
be varied by 16 stepping jacks and each wall is equipped vfith 58 central 
pressure pickups (dia. = 0,4 mm) and a few lateral pickups. The average 



test conditions are: PT = 1.7 bars; TT = 290°K; RC = 10^ (c = 400 mm; Mo = 
0.84) . 

The C5 body is an axisymmetric model with a slowly variable section 
approximating the outside dimensions of a transport aircraft around its 
centerline. The model used was 400 mm long (fig. 1). It is made of 
3-mm-thick steel parts generated by rotation, which, if necessary, could 
enable rapid cooling during a cryogenic gust. Pressure distribution is 
recorded by 45 0 . 3-mm-dia . pickups aligned along a generatrix. In 
addition, with x/c = 2.1%, 4 pickups are arranged in a cross check for proper 
alignment of the model in relation to the tunnel axis. The model is held 
at the rear by a freely rotating rod (fig. 1). A cable system driven by a 
stepping motor turns the model about its axis and therefore the line of 
pressure pickups as well. The wind tunnel control computer effects several 
successive rotations of 45° or 90° separated by pressure data acquisition. 
Further to the downstream side a set of 2 swivel joints enable the C5 body 
to be set at an incidence angle at any height in the tunnel. Boundary 
layer transition is triggered by a plastic strip approximately 0.2 mm thick 
glued at x/c = 7%. 
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3 - ADAPTATION METHOD 

The physical shifting of the top and bottom flexible walls is 
two-dimensional (ref. 6). But the three-dimensional nature of the flow is 
tqken into account by the adaptation calculation program. 



3.1 - Principle 


The adaptation calculation is based on a schematic representation of 
the model by singularities and on the contribution of the images of these 

singularities in relation to the walls of the channel. 

t 

The model is represented schematically by a distribution of 3D sources 
and horse-shoe vortexes (of infinitely small width) placed on the axis of 
the model (fig. 2). The intensities of these singularities are deduced 
from the pressure measurements made on the top and bottom walls. 

The interference of the walls, represented by the disturbances due to 
the images alone, is then calculated at the model (fig. 2). This 
interference is finally canceled by superimposing an opposite interference 
due to an estimation of a new shape of the walls. This new shape is a 
combination of a symmetrical deformation (due to the sources) and an 
antisyrnmetrical deformation (due to the vortexes). 

In reality this process is simplified so that we can go directly from 
the peripheral pressure measurements t.o the new shapes thanks to two linear 
operators connected to the sources and the vortexes respectively, which 
depend on the geometry of the tunnel (fig. 2). 

3.2 ~ Possibilities 

The use of image singularities is not easy around flexible walls of any 
geometric shape. Consequently, the first version of this method required 
an initial tunnel with flat parallel walls and calculating of a "final 
shape" in one step. 

The present version enables us to start with any shape of walls. An 
initihi treatment at the peripheral pressure measurements replaces the real 
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tunnel with a parallelepiped control surface. The normal process is then 
applied . 

Under these conditions, since the adaptation is theoretically not 
achieved after one iteration, it is possible to follow with several 

iterations taking into account the successive shapes assumed by the walls. 
In this way the residual interferences are minimized much better than with 
the first version by coming closer to the perfect shapes adapted to the 
linearity hypotheses in the calculation, to the errors due to 
discreteness, and to the measurements and mechanical precision. Furthermore, 
the Mach number Mo measured to the upstream side of the channel and the 
true Mach number of the test related to the flow around the model are 
closer and closer after each repetition. 
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The advantage of this second version of the adaptation method enables a 
test to be carried out starting with a partially unblocked tunnel. This 
configuration becomes necessary at high Mo Mach values together with a 
model of large dimensions. 

3.3 - Remarks 

In all cases, with or without adaptation, the shapes of the top and 
bottom walls take into consideration on a global basis the viscous effects 
of the boundary layers which develop on the four sides of the tunnel. The 
shape of the top wall, for example, takes into consideration the boundary 
layer which develops on it (calculation from the pressure measurements) and 
the boundary layer on one of the lateral walls (calculated as a flat plate 
at Mo)i To do this, the sum of the two traveling thicknesses of these 
boundary layers is added to the perfect fluid boundary coming from the 



adaptation calculation. 


The compressibility effect is taken into account by the coefficient 
] 

which comes into play in the Prandtl-Glauert rule. 

4 - RESULTS 


The results can be subdivided into three parts: 

t 

- Zero incidence 0.6 <. Mo _< 0.95 

~ Incidence 0( ; = 8°; Mo = 0.84 

- Zero incidence, model off center. 

4.1 - Adaptation at Q{. - 0°, 0.6 < Mo £ 0.95 

Mo = 0.6 At Mo = 0.6 the mark of the model on the unadapted walls Is 
quite weak (fig. 3). The swivel joint on the downstream side creates a 
speed disturbance of the same order of magnitude. Adaptation is 
performed in a single iteration. We can observe (fig. 3) that the 
adaptation around the swivel joint located at the ends of the deformable 
walls is approximately correct (model and swivel joint dimensions close, 
Mach numbers similar). On the model (fig. 4) the Mach number distribution 
undergoes an overall decrease which is greater in the high speed zones. 
Figure 5 indicates marked unblocking of the walls. 

Mo “ 0.7 Adaptation at Mo = 0.7 (fig. 6) constitutes an initial point 
of comparison. Figure 7 shows the sharp correction in the distribution of 
Mach number obtained in one iteration. Figure 8 indicates good agreement 
between the result of the adaptation performed at T2 , that of the strictly 
3D adaptation performed in Berlin (TUB ref. 7), and the measurements from a 
1 a rg e Wind tunnel (NASA-Atnes braking coefficient = 1.4%). 



Mo ~ 0.75 Mo = 0.75 seems to be the limit of adaptation in a single 


operation for the configuration considered (figs. 10 and 11). 
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The walls are adapted as of the first iteration around the model, but a 
second iteration fine tunes the adaptation around the swivel (fig. 12). 

Mo ~ 0.84 At Mo - 0.84 the mark of the model between non-adapted walls 
starts to become large (fig. 13). Further to the downstream side the 
influence of the swivel is marked by a very strong pressure disturbance 
which induces bulging on the upstream part of the walls. One also observes 
a kind of positive average speed gradient (fig. 13) which disappears upon 
adaptation . 

After, for Mo > 0.84, the distribution of the Mach number on the model 
is plotted with reduced variable M/Mo to bring down the deviations due to 
slightly different values of Mo. 

In figure 14 we see the sharp decrease in the Mach number on the model 
between the non-adapted case and the first iteration. This tendency is 
mostly noticeable at the maximum cross-section where the impact intensity 
greatly decreases, and on the rear portion of the model where the viscous 
effects must change. 

We find it necessary to perform three iterations to converge toward a 
stable solution. Nevertheless, the influence pf the first iteration is by 
far the most important (figs. 14 and 15). 

A comparison with the results of a 3D adaptation conducted in Berlin 
(TUB ref. 7) shows the effectiveness of this method. .. 
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Mo i=s 0.9 In order to avoid blocking of the channel in the vicinity of 
the model, this test starts off using a shape more or less adapted at Mo - 
0.89 as the initial shape. Under these favorable conditions, a single 
iteration is sufficient to adapt the tunnel (fig. 16), as demonstrated by 
the identical results of iterations 2 and 3. 

Comparison with the results from NASA-Ames indicates a very good 
agreement -even at the impact points- very noticeable at Mo = 0.9 (fig. 

170 . 

Mo 1=1 0.95 At Mo ~ 0.95 we reach the usage limit of the present method 
as applied to wind tunnel T2 and the particular model. 

The speed disturbances on the top and bottom walls at the maximum 
cross-section of the model become substantial and the Mach number is 
greater than 1 (fig. 18). The connection between these disturbances and 
the calculated shape of the waifs intensifies and convergence of the 
adaptation process is more and mpre difficult to achieve. 

In addition, during this test series, the absence of automatic Mach 
number regulation of the test caused some slight variations in this 
parameter from one gust to another, that is, from one iteration to the next 
(delta Mo < 0.01). These changes in Mo complicate the relationship between 
disturbance and wall shape and increase the difficulty of convergence. 

At about Mo ~ 0.94, three successive iterations seem to approach a more 
or less converged solution; but the variations in Mo from one iteration to 
another (0.935 _< Mo < 0.944) prevent proper stabilizing of the process. 
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When trying to reach Mo = 0.95, the disturbance/wall shape relationship 

JL 

is so great that the adaptation process no longer converges. , This results 
in strong disturbance zones along the walls which vary in position and 
modulus. On the model the longitudinal movement of the impact around the 

maximum cross-section is great, and a very unstable supersonic pocket forms 
behind it. 

Nevertheless,' a comparison was made at Mo = 0.95. For this purpose one 
iteration was chosen among the five successively run iterations. The 
selection criterion was based on seeking the best coherence between the 
shapes of the walls at the various Mach number values Mo. This comparison 
with the results of the Berlin TU and the NASA Ames indicates good overall 
agreement on the model (fig. 19). Figure 18 shows the strong speed 
disturbances on the walls around the model and the swivel. 

Ob s ervati on : When Mo rises, the linearity hypothesis contained in the 
adaptation calculation grows less and less justified. Nevertheless, good 
coinciding at Mo = 0.90 (fig. 17) would indicate that this hypothesis only 
goes far out of line when very close to Mo - 1. 

4.2 - Adaptation at &(- 8°. Mo - 0.84 

This configuration comprises three tests: 

- Measurements between diverging walls compensating 
only for the boundary layers (non-adapted case). 

- First iteration of adaptation, 

! 

- Second iteration of adaptation. 
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In these three tests, the C5 model is shifted in incidence. The degree 
of incidence, recorded with a cathetometer in the presence of the gust, is 
about 8°, give or take 5'. Actually, in order to prevent the 

combination of an incidence problem with an off-center problem, the model 
is inclined by rotating it about its center. 
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Figure 20 shows the small change made in the disturbances on the 
non-adapted walls by giving the model and angle of incidence. Let us point 
out the presence of a longitudinal pressure gradient in the tunnel, as in 
the previously planned non-adapted cases. The Mach number of this test is 
taken to be equal to 0.817, but its value at the model is closer to 0.84. 

Furthermore, during each of these three gusts, the C5 model is rotated 
about its longitudinal axis in order to make measurements every 45°. It 
seems that this rotation is affected by an initial error: the first 45° 
rotation is only partially performed (about, halfway) perhaps because a 
force is applied to the driving cable without there occurring any actual 
rotation. This phenomenon varies slightly from one test to another. It 
can be brought to light by plotting the nine dis tribution^ of reduced Mach 
nurnbeb Of the non-adapted case (fig. 21); this initial deviation is clear 
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in the central zone (x/c = 50%) where the shapes of the curves vary 
tremendously. A similar deviation, but. probably of lesser magnitude, is 

I 

found in the adapted case (first iteration, fig. 22). 

On figures 21 and 22, let us point out the large variation in the Mach 
number distribution depending on the generatrix. Near the maximum 
cross-section we observe a noticeable rise in the impact and a decrease in 
its intensity when the generatrix goes from the "lower surface" to the "top 
skin." Likewise, shifting of the stagnation point results in a 
considerable dissymmetry in the flow depending on the angle, near the nose 
of the model (0% < x/c £ 15%). 

Figure 23 shows another representation of a slice of depending on the 
polar angle. This is a plot of the Mach number measured by the four 

i 

pickups arranged in a cross in the section x/c = 21%. When the model 
rotates about itself the measurements of the four pickups change in harmony 
depending on the angle; these curves indicate relatively good right/left, 
symmetry in relation to a vertical axis corresponding to Saxe > 180 °* The 
less sophisticated finish of the three pickups added in this section (£7, 

X 5 $) could explain the differences in levels between the curves. 

As concerns the adaptation itself, figure 2 A shows a change in the Mach 
number on the inclined model. It was found necessary to perform two 
iterations to obtain a stable result. An overall decrease characterizes 
the adaptation of the top and bottom walls, but no other result coming from 
another facility enables the quantitative validity of this converged 
configuration to be verified. 



5 - MODEL OFF CENTER 


Since the new adaptation calculation theoretically enables the model to 
be set off center in the tunnel, we tried to test this possibility. In 
order to do so, the model was placed level in the middle of the tunnel at 
-57% with zero incidence. 



Key to figure 

1: centered 2: off center 

As of the first iteration calculation of the linear operator /}( fig . 2) 
is defective. Remember that this operator, pertinent to the schematic 
representation according to sources (symmetrical part), enables direct 
transition from the pressure measurements on the walls to their symmetric 
strain. In the present case, the discrete solution JjL(x) is found for an 
approximately correct average behavior, but it is modulated by a 
substantial oscillation which forbids its use. Indeed, the multiplication 
of these oscillations by the pressure disturbances measured on the walls 
causes breaks in the calculated shapes of the walls. Wham's more, the 
average curve is too imprecise to be substituted for the discrete solution. 



This problem seems to be purely mathematical and further study of this test 
situation is planned. 

The of f -centering imposed here, however, is very important. We checked 
by calculation that a smaller off-centering ~ 10% leads to a correct 

mathematical solution of the operator , and therefore for the wail 
shapes. This configuration has not been subjected to experimentation. 

6 - CONCLUSIONS 


These tests enabled to validate the implementation of an iterative 
adaptation process around a large axisymmetric model at the T2 wind tunnel. 


The adaptation calculation is based on a' calculation of wall 
interferences in three dimensions. The advantage of it is that it enables 
any initial shape of the flexible walls to be taken into account, and 
provides for adaptation around an axisymmetric model at an angle of 

i 

incidence. It can theoretically deal with cases of models off center in 
the tunnel. 

Movement of the flexible walls is still two-dimensional. 


14 

With zero incidence and the model centered one adaptation iteration 
seems to be sufficient for Mo < 0,84. Beyond, two to three iterations' are 
necessary, but the first one provides for most of the correction. In these 
iterative cases, convergence of the process is effective up to Mo =■ 0.9 
included; for Mo > 0.94 the shapes of the walls oscillate without 
converging . 

! 
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Comparisons with an entirely three-dimensional adaptation conducted in 
Berlin (TUB) and with tests conducted in a large tunnel (NASA Ames, 11 ft. 
x 11 ft., braking coefficient = 1.4%) reveal good effectiveness at the 
model level (braking coefficient 1.8%) of the pseudo three-dimensional 
adaptation of the T2 wind tunnel; the residual interferences are very 
small . 


At Mo = 0.84 and < A = 8° of incidence, the adaptation process converges 
after three iterations towards a likely solution, 


A single of f -centering test of the model, which was probably too great 
(fi/T* did n °t enable validation of the method for this type of 

configuration. Mathematical calculation of the linear symmetric operator 
(see § 3.1) gives oscillating solutions which cannot be used for proper 
determination of the wall shapes- 
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Figure 1: T2 wind tunnel 
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Key to figure 1 

1: compressor 2: reservoir: 80 bars 3: heater 4: nitrogen 5: 
adaptable walls 6: C5 body, C: 400 mm, max :58mm 7: Incidence control 
joints 8: transverse rotation cable 
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"3D" Adaptat 


ion at tunnel center - Figure 2 
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Key to figure 2 

1: 2D tunnel 2: walls 3: variable top and bottom 4: stationary sides 
5: movement of top and bottom ~ symmetric movement + antisymmetric 
movement 6: model represented by [vortex sources] 7: 3D on the tunnel 
axis 8: measured 9: wall interaction 10: model 11: Method (example 
sources = > symmetric movement) 12: induced 13: geometry 14: linear 
operator 15: with 16: To cancel the interaction of the walls u^(x) on 
the tunnel axis, -ui(x) is created by symmetrical movement of the walls 
17: Procedure similar for antisymmetric movement with vortexes and 



3 - Mach number on the flexible walls 
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Key to figure 3 

1: top wall 2: bottom wall 3; non-adapted 4: adapted 
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Figure 5 - Shape of flexible walls 



Key to figure 5 
1: non-adapted 2: adapted 
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Figure 6 - Mach number on the flexible walls 
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Key to figure 
1: top wall 


6 

2: bottom wall 


3: non-adapted 


4: adapted 
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Figure 7 - Mach number on the C5 body 



Key to figure 7 
1: non- adapted 


2 : adapted 3: adapted 3D-T 
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Figure 8 - Mach number on the C5 body 



Key to figure 8 

1: adapted 2: 3D adapted 3: vertical generatrix 4: horizontal 
generatrix 
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Figure 9 - Shapes of flexible walls 



Key to figure 9 
1; non-adapted 2: adapted 




Figure 10 ~ Mach number on the flexible walls 



Key to figure 10 

1: top wall 2: bottom wall 3: non-adapted 4: adapted 
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Figure 12 ~ Shapes of flexible walls 
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Figure 13 - Mach number on the flexible walls 



Key to figure 13 

1: top wall 2: bottom wall 3: non-adapted 4: adapted 
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Key to figure 14 

1:^ non-adapted 2: first iteration 3: third iteration 4: 3D adapted 


oo 





BOO 


Figure 15 - Shapes of flexible walls 



Key to figure 15 

1: non-adapted 2: first iteration 3: third iteration 
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Key to figure 16 

1: top wall 2: bottom wall 3; adapted 
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Figure 17 -- Mach number on the C5 body j 



Key to figure 17 

1: adapted 2: vertical generatrix 3: horizontal generatrix 
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Figure 18 -• Mach number of flexible walls 



Key to figure 18 

1: top wall 2: bottom wall 
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Figure 19 - Mach number on the C5 body 
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Key to figure 19 

1: adapted 2: 3D adapted 

generatrix 
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3: vertical generatrix 4: horizontal 
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Figure 23 


- Rotation of 4 pickups of section x/c-217o 





Key to figure 23 

1: rotation positions 2: four pickups of the section x/c=21% 
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Figure 25 - Adaptation - Shapes of flexible walls ( =8°) 
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key to figure 25 

1. non- adap t.ed 2 : first iteration 3; second iteration 
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resume O'auteur : Lors de ces essais, Tes parois bidimensionnelles haute et basse de lai 

soufflerie 72 sont adaptees iterativement. Le calcul d'adaptation tient compte de la 
tridimensionnal its de 1 'ecoulement ; 4 chaque itEration, il reduit 1 1 interference des 
quatre parois au niveau de la maquette. Cette mSthode permet done de partir des formes 
de parois quelconques ; un second avantage consiste 4 adapter les parois autour d'une 
maquette portante. 

Les essais ont et§ realises avec une maquette axisymetrique C5 
(biocage 1,8 X L ongueur 400 mm) , 4 temperature ambiante et pression d'arret moyenne de 
1,7 bar (Nombre de REYNOLDS base sur la longueur du corps : 10 7 4 Mo = 0,34) . Ilscompren- 
nent : - un balayage en nombre de Mach entre 0,6 et 0,95 4 incidence nulle 

' - une adaptation 4 incidence o< * 8° 4 Mo ■ 0,84 
- un essai non converge avec le modeie fortement decentre, 4 

incidence nulle. 

Pour certaines valeurs de Mo, 4 incidence nulle, des comparaisons 
sont faites avec desresultats d‘une vraie adaptation tridimensionnelle (TU BERLIN) 
et de tests effectuSs dans une grande soufflerie (NASA Ames). 


NOTIONS O'lNOEXAGE : 

Corps C5 - Adaptation de parois - Adaptation en transsonique - Interference 
tridimensionnelle. 
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NOTATIONS 


B Largeur 


veine “ 390 nan 


C Longueur modele “ 900 mm 

H Hauceur veiue * 370 mm 


K P (x_})= 4 Tr m n m* B*+(n-^f H 4 ]]* 


Coefficient d' influence dcs 
singular ites au niveau de3 
prises de pression parols 
flexibles. 


K i (x -^ )= 4TT [(x-jff^(m‘B 1, +’n‘H*)] vl 

(m,n) *(0,0 ) P 


Coefficient d ' influence des 
singularity images seules 
au niveau de l'axe veine. 


Nombre de Mach infini de 1'essai 
Pression totale 

Nombre de REYNOLDS de l'ecoulement (lie & C) 


Temperature totale 

Coordonnees cartesiennes 
repSre orthonorme 


2 D Abreviatiou de bidimensionnel 

3 D Abreviatiou de tridimensionnel 

0( Angle d' incidence affiche avant essai 

(9 - \J 1 - Mo 2 , coefficient de compressibilite 


gene, verti. 
gene, horiz. 



Essais realises 3. la Technische Universit^t 3erlin 

Generatrice verticale ) Eg3al3 >, ASA ^es 
Generatrice horizontale ) 
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fluence des 
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I - INTRODUCTION 


Une premiiire campagne d'essais (septembre 1984, ref. / 1 / ) avail: 
montre la validita d'ane adaptation bidimensionnelle auCour de modules 
tridimensionnels. La methode d'adaptation s'appuyait sur une correction de 
vlcesse calculee an un pas, en tridimensionnel. 

La presente campagne prolonge la precedents et etend le domaine de 
validite du processus employe grace aux caractdri3tiques suivantes : 

- L'obstacle considSre est le corps axisymetrique C5 (ref. / 2/), 
dej4 teste en 1984. Mais la maquette est, cette fois, plus volumineuse (C • 400 mm ; 
coefficient de blocage 1,8 Z) . 

- La methode de calcul, mise au point par L. LAMARCHE (ref. /3/) 
differe de la premiere version utilisae : il s'agit d'une methode iterative 
plus precise et qui permet de corriger les interferences de parois autour de 
modSle en incidence ou de module decentre. 


2 - SOUFFLERIE T2 ET MAQUETTE 05 


La soufflerie T2 est une installation cryogenique 4 circuit ferme, 
fonctionnant 1 induction, par rafales de 1 4 2 minutes (fig. 1) - Ref. / 4/ et 
/5/. L'actuelle campagne d'essais a ete rdalisee uniquement 4 temperature ambiante. 

La veine d'essai (H ■ 370 mm, B » 390 mm 4 1' entree) comporte deux 
parois flexlbles haute et basse ; chacune de ces parois esc deformee par 16 
verins pas 4 pas et est equipee de 58 prises de pression centrales (0 ” 0,4 mm) 
et quelques prises laterales. Las conditions moyennes d'essai sont : PT ”1,7 bar ; 
TT- 290 °K ; RC » 10 7 (c - 400 mm ; Mo - 0.34). 

Le corps 05 est un modSle axisymetrique de section lentement variable, 
schematisant 1 ' enccmbremenc d'un avion de transport ramene autour de son axe. La 
maquette utilisee a 400 mm de long (fig.. 1). Elle est constltuee de pieces de 
revolution en acier de 3 mm d'epaisseuc, permectant eventuellement un refroi- 
dissement rapide en cours de rafale cryogenique. La distribution de pression est 
relevee par 45 prises da 0 ■ 0,3 mm, alignees sur une generatrlce. De plus, en 
x/c « 21 Z, 4 prises disposees en croix sont destinees 4 verifier le bon alignement 
de la maquette suivant l'axe veine. Ce mod41e esc tenu par l'arriere par une 
canne (fig. 1) libre en rotation. Un systSme de cable entraine par un moteur pas 

4 pas permet alnsi de falre tourner la maquette autour de son axe et par suite 
la gSneratrlce de prises de pression. L'ordinateur de commande de la soufflerie 
impose plusieurs rotations successives de 45° ou 90°, separees par une acqui- 
sition de pression. Plus en aval, un jeu de 2 rotules permet de placer le corps 

05 en incidence 4 n'importe quelle hauteur dans la veine. La transition de la 
couche limite esc declenchee par une bande plastique d'environ 0,2 mm d'epaisseur, 
collee en x/c * 7 %. 
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3 - METHODE D' ADAPTATION 


La deplacement physique des parois flexibles, haute ec basse, eat 
bidimensionnel (Rdf. /6/). Mala le caractere tr idimensionnel de l'ecouiement 
est pria en compta par le programme de calcul d'adaptation. 

3.1 - Principe 

Ce calcul d'adaptation est base sur une schematisation de la maquette 
par des singularites et 3ur la contribution des images de ces singularites par 
rapport aux parols de la veine. 

La maquette est schematises par une distribution de sources 3D et de 
tourbillons en fer 3 cheval (de largeur infiniment pecite) placee sur l'axe du 
models (fig. 2). Les intensites de ces singularites sont deduites des mesures 
de presslon faltes sur le3 parols haute et basse. 

L' interference des paroi3, representee par les perturbations dues aux 
images seules, est ensuite calculee au niveau du models (fig. 2). Cette inter- 
ference est enfin annulee par une superposition d'une interference opposee, due 
3. une estimation d'une nouvelle forma des parois. Cette nouvelle forme est une 
combinaison d'une deformation symetrique (due aux sources) et d'une deformation 
antisymetrique (due aux tourbillons). 

En realite, ce processus est simplif lee de maniSre 1 passer directement 
des mesures de presslon parietale aux nbuvelles firmes, grace 3 2 operateurs 
lineaires respect tveaent lies aux sources et aux tourbillons et dependant de la 
geometrie de la veine (fig. 2). 


3.2 - Possibilites 


L'utilisation des singularites images n'est pas facile autour 
de parois flexibles de geometrie quelconoue. Aussi, la premiere version de 
cette methods exigeait une veine de depart 3 parois planes paralleles et calculait 
une "forme finale" en un pas. 

La version prSsente permet de demarrer de formes de parols quelconques. 

Un premier traitement au niveau des mesures de presslon parietale remplace la 
veine reelle par une surface de controle parallelipipedique. Ensuite, le 
processus normal est applique. 

Dans ces conditions, 1' adaptation n'etant theoriquement pas realisee 
aprSs une iteration, il est possible d'enchainer plusieurs iterations qui 
tiennent compte des formes successlves prises par les parois. On minimise alnsi 
beaucoup mieux qu'avec la premiere version les interferences residuelles en 
3'approchant de plus pr3s des formes adaptees parfaites aux hypotheses de 
linearlte dans le calcul, aux erreur3 dues 3 la discretisation et aux mesures 
et 3 la precision mecanique prSs. D' autre part, le nombre de Mach Mo mesure 
3 l'amont de la veine et le vrai nombre de Mach de l'essai lie a l'ecouiement 
autour de la maquette sont de plus en plus voisins au fur et 3 mesure des itera- 
tions. 
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L'avantage de cette deuxi&ne version de la methode d 'adaptation permet 
egalement d'effeccuer un essai en partant d'une veine partiellement debloquee. 
Cette configuration devient necessaire 3 grand nombre de Mach Mo en presence 
d’une maquette de fort encombrement. 


3.3 - Remarques 


Dans tous les cas, adaptes ou non, les formes des parois haute et 
basse tiennent globaleaent compte des effets visqueux des couches limites qui 
se developpent sur les quaere cotes de la veine. La forme de la paroi haute, 
par exemple, prend en consideration la couche limite qui se developpe sur elle 
(calcul 3. partir de3 mesures de pression) et la couche limite exi3tant sur une 
des parois lacera.les (calculee comme une plaque plane 3 Mo). Pour cela, on 
rajoute 3 la frontidre fluide parfait issue du calcul d'adaptation la somme des 
deux epaisseurs de deplacement de ce3 couches limites. 

L'effet de compressibilite est pris en compte par l'intermediaire 
du coefficient (3 '/ 1 — M q intervenant dans la rlgle de PRANDTL - GLAUERT . 


4 - RESULTATS 


L'ensemble des resultats peut etre subdivise en trois parties : 

- Incidence nulle 0,6^ Mo ^0,95 

- Incidence » 8 3 ; Mo « 0,84 

- Incidence nulle, modSle decentre. 


4.1 - Adaptation 3 o< » 0*. 0,6^ Mo C 0,95 

Mo » 0 , 6 1 A Mo »0,6 , la trace de la maquette sur les parois non 

adaptees esc ass ez faible (fig. 3). La rotule, en aval, cree une perturbation 
de vltesse du meme ordre de grandeur. L'adaptatlon esc reallsee en une seule 
iteration. On peut remarquer (fig. 3) que 1' adaptation, aucour de.la rotule situee 
3 1'extremitS des parois dSformables, est relativement correcte (encombrements 

modele - rotule voisins, niveau de nombre de Mach voisins). Sur la maquette 

(fig. 4), la distribution du nombre de Mach subit une decroissance d'ensemble, 
plus importance dans les zones de vitesse elevee. La figure 5 indique un net 
deblocage des parois. 

Mo ■ 0,^1 L' adaptation 3 Mo “0,7 (fig. 6) const itue un premier 
point de comparaison. La figure 7 montre la necte correction de la repartition 
du nombre de Mach obtenue en une Iteration. La figure 8 indique le bon accord 
entre le resultat de l'adaptation effectuee 3 T2, celui de l'adaptation rigou- 
reusement 3D realisee 3 Berlin (TUB ref. / 7 / ) et les mesures relatives 3 une 
grande soufflerie (NASA-Ames coef. blocage * 1,4 7 ..). 

Mo “ 0 , 7 5 | Mo “ 0,75 semble constituer la limite de l'adaptation en 
une seule iteration, dans la conf iguraiton consider.ee (fig. 10 et 11). Les 
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parols sonc adaptSes d§s la premiSre iteration aucour da la maquette, mala 
uae deuxieme iteration affine l'adaptatlon autour de la rotule (fig. 12). 

• lte'^0^84^ A ® « 0,84, La trace de la maquette entre parols 

non adaptees commence 4 etre importance (fig. 13). Plus en aval, l'influence 
de la rotule est marqude par une trSs forte perturbation de pression, qui 
induit des bombements sur la par tie amonc des parols. On remarque egalemenc une 
sorte de gradient moyen positif de vitesse (fig. 13), qui diparalt 4 l'adaptatlon. 

Par la suite, pour Mo ^0,84, la repartition du nombre de Mach, sur 
la maquette sera traces en variable reduite M/Mo pour diminuer les icarts dus 
a des valeurs legSrement differences de Mo. 

Sur la figure 14, on observe la nette decroissance du nombre de Mach 
sur le module, entre le cas non adapte et la premiere Iteration. Cette tendance 
est surtout sensible au maltre couple, ou l'intensite du choc decroit fortement, 
et sur la partie arriSre de la maquette, oil les effets v.isqueux doivent evoluer. 

On. constate la necessite d'effectuer trois iterations pour converger 
vers une solution stable. Toutefois, l'influence de la premiSre iteration est 
nettement la plus importance- (fig. 14 et 15). 

Une comparaison avec les resultats d'une adaptation 3D realises 4 
Berlin (TUB ref /7/)'moncre l'efficacite de la methode presente. 


j Mo » 0 , 9 1 Pour eviter un blocage du canal au niveau du modSle, cet 
essai demarre en uti’lisant une forme a peu prSs adaptee 4 Mo » 0,89 comme forme 
initiale. Dans ces conditions favorables, une seule iteration suffit 4 adapter 
la veine (fig. 16), comme l'ont raontrd les resultats identiques des iterations 
2 et 3. 


La comparaison avec les resultats issus de la NASA - Ames indique un 
accord trSs correct, meme aux niveau du choc, cres sensible 4 Mo « 0,9 (fig. 17). 

Mo » o', ""sj A Mo * 0,95, nous acteignons la limite d'utilisation 
de la presence metKode appliquee 4 la soufflerie T2, et 4 la maquette consi- 
deree. 



Les perturbations de vitesse sur les parois haute et basse, au droit 
du maltre couple du modSle, deviennent importances et le nombre. de Mach est 
superieur 4 1 (fig. 18). Le couplage entre ces perturbations et la forme calculee 
des parois devient plus fort et la convergence du processus d 'adaptation est 
de plus en plus difficile 4 realiser. 

D'autre part, durant cette campagne, l'absence de regulation automa- 
tique du nombre da Mach de 1' essai entralne quelques legSres evolutions de ce 
paramStre d'une rafale 4 l'autre, c'est-4-dire d'une iteration 4 la suivante. 
(AMo^ 0,01). Ces variations de Mo compliquenc le couplage perturbation -forme 
de parois et augmentent la difficulte de convergence.. 

Autour de Mb- ■ 0,94, trois iterations successives semblent cendre 
vers une solution plus ou molts convergee ; mais les variations de Mo d'une 
iteration 4 l'autre (0, 935<Mo ^ 0 ,944) empechent une stabilisation correcte 
du processus. 
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Lorsqu'on chercha 1 atteindre Mo ■ 0,95, le couplage /perturbation - 
forme parois/ est si fort que le processus d'adaptation ne converge plus. 

Cela se traduit au niveau des parois par des zones de force perturbation qui 
varient en position at en mbdule. Sur la maquette, le deplacement longitudinal 
du choc autour du maltre couple est important et derriire lui se reforme une 
poche supersonique cres instable. 

Une comparaison a tout de meme ete faite 1 Mo * 0,95. Pour cela, 
une iteration a eta choisie parmi cinq iterations successivement soufflees. Le 
critere de ce chotx est fonde sur la recherche de la meilleure coherence entre 
les formes de parois, aux divers nombre de Mach Mo. Cette comparaison avec les 
resultats du TU Berlin et de la NASA Ames indique une bonne concordance 
d 'ensemble sur le module (fig. 19). La figure 18 montr’e les fortes perturbations 
de vitesse sur les parois, au niveau du modele et de la rotule. 

Remarque : Lorsque Mo crolt, l’hypoth&se de linearite contenue dans 
le calcul de 1 'adaptation devient de tnoins en moin3 Justifiee. Toutefois la 
bon recoupemenc 1 Mo * 0,90 (fig. 17) semble indiquer que cette hypothese ne 
devient vraiment abusive que trSs pres de Mo “ 1 . 

4.2 - Adaptation a C< ■ 8 a . Mo » 0,84 

Cette configuration comprend trois essais : 

- Mesures entre parois divergences compensant juste les couches 
limites (cas non adapte) . 

- lere iteration d 'adaptation. 

- 2Sme iteration d 'adaptation. 

Dans ces troi3 essais, le module C5 est mis en incidence. La valeur de 
1' incidence, relevee au cathetomStre pendant la rafale, est d' environ Of » 8°- 
a 5' pres. En fait, pour ne pas meler un propleme d' incidence et un problSme 
de decentrage, la maquette est mise en Incidence par rotation autour de son 
centre. 
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La figure 2Q montre la faible changement apporte au niveau des 
perturbations sur las parois non adaptSes par la tnise en incidence de la 
maquette. Remarquons la presence dans la veins d'un gradient longitudinal de 
pression, comae dans les cas non adaptes precedemment envisages. Le nombre de 
de Mach de cet esaai. est pris egal a 0,817 mais sa valeur au niveau du modele 
est plus proche de 0,84. 

D'autre part, durant chacune de ces trois rafales, le modSle C5 est 
mis en rotation autour de son axe longitudinal afin d'effectuer des mesures 
tous les 45*. II aemble que cette rotation soit entachee d'une erreur initiale : 
le premier pas de 45* n'est executS qu'en partie (environ la moitiS), peuc-etre 
I. cause d'une misci en tension du cable d'entrainement sans rotation effective ; 
ce phenomene varia leg&rement d'un essai a l'autre. On peut le mettre en 
evidence en traqant les neuf distributions du nombre de Mach reduit du cas 
non adapte (fig. 21) ; ce decalage initial est nec dans la zone centrale (x/c ■ 50 2) 
oil la forme des courbes evolue enormement. On retrouve un decalage analogue, mais 
vraisemblablement moins important dans le cas adapte (premiere Iteration, fig. 22). 

Notons, sur ces figures 21 et 22, la grande evolution de la distri- 
bution du nombre de Mach en fonction de la generatrice. PrSs du maltre couple, 
on' observe une avane.e sensible du choc et une decroissance de son intensite 
lorsque la gSneraerice passe de "1 ' intrados k l'extrados". De meme, le deplacement 
du point d'arret se traduit par une dissymetrie notable de l'ecoulement en 
fonction de 1'ang.le, pr3s du nez du modSle ( 0 2^x/cj^15 2). 

La figure 23 montre une autre representation d'une tranche d'Scou- 
lement en fonction del'angle polaire. II s'agit du trace du nombre de Mach 
mesure par les quatre prises placeesen croix dan3 la section x/c » 21 2. Lorsque 
le modele tourne sur lui-meme, les mesures des quatre prises evoluent de faqon 
analogue en fonction de l'angle ; ces courbes indiquent une symetrie gauche - 
droite relativement correcte, par rapport a un axe vertical correspondant & 

0axe > 180*. L'Stat de finition moins elabore des trois prises rajoutees dans 
cette section ( 0 , x , * ) peut expliquer les decalages de niveau existant 
entre les courbes. 

En ce qui concerne 1' adaptation elle-meme, la figure 24 montre une 
evolution du nombre de Mach sur la maquette en incidence. On constate la necessite 
d'effectuer deux itSrations pour obtenir un resultat stable. One decroissance 
d'ensemble caracterise l'adaptation des parois haute et basse, mais aucun autre 
resultat issu d'une installation diffSrente ne permet de controler la validite 
quantitative de cette configuration convergee. 


5 - MODELE DECENTRE 


Le nouveau calcul d' adaptation permettant en theorie de decentrer le 
modlle dans la veine, nous avons essaye de tester cette possibilite. Pour cela, 

la maquette a ete placee dans la plan horizontal median, a “ -57 Z, k 

incidence nulle. 
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DSs la pre.mitre iteration, le calcul de l'optrateur lintaire Q (fig. 2) 
est en defaut. Rappelons qua cet operateur, relatif i la schtmatisatioh par 
sources (oartie symetrique) permet de passer directement dea mesures de pression 
sur les parois a la dtformte symetrique de celles-ci. Dans le cas present, la 
solution discrete Q (x) trouvee a une allure moyenne approximativement correcte 
mai3 elle esc modules par une oscillation importance que interdit son utilisation. 
En effet, la multiplication de ces oscillations par les perturbations de pression 
mesurees sur les parols provoque des cassures sur les formes de parois calcultes. 
D'autre part, la courbe moyenne est trop imprecise pour etre substitute I la • 
solution discrete. Ce probltme semble etre purement mathematique et une etude 
plus approfondie de ce cas d'essai est envisagee. 

Toutefois, le decentrement impost ici est tres important. Sous avons 

verifit par le calcul qu'un dtcentrement plus faible ■ - 10 % conduit 

a une solution mathtmatique correcte pour l'optrateur Q , done pour les formes 
de parois. Cette configuration n'a pas ttt exptrimentalement rtalisee. 


i une 

i ntcessite 

ssance 6 - CONCLUSIONS 

m autre 

ilidite 

Cet ensemble d'essais a permis de valider la mise en oeuvre, i la 
soufflerie T2, d'un processus ittratif d'adaptation autour d'un modtle axisymt- 
crique de grande taille. 

Le calcul d'adaptation. est fondt 3ur un calcul d' interferences des 
parois en tridimensionnel. II a pour avantages de pouvoir prendre en compte 
des formes initiales quelconques des paroi3 flexibles, et d'etre capable d' adapter 
'er le autour d'un modtle axisymetrique en incidence. II peut thtorlquement traiter 

• cela, j des cas de maquettes dtcentrees dans la veine. 

a Le dtplacement des parois flexibles reste bidimensionnel. 
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A incidence nulle, module centrS, une iteration d'adaptation semble 
! ✓ n 84 Au-dell, deux 1 trois iterations sont necessaires, 

'“I. U P««» d. U correction. Dana ... ... 

ltfiratlfs la conversance du proceaaua aac effective jnaqn 4 Mo - 0,9 conprla , 
pour Mo >0,94, les formes des parols oscillent sans converger. 

Des comparaisons avec une adaptation complement tridimenslonnelle 
r dal la fie > Berlin (TO) « a~ “US Xlu, 

r“v^*d. 1. aecnect. (bloc,,. 1.8 p.jndo rrldMenalonnell. 

de la soufflerie X2 ; les interferences residuelles sont tres faibles. 

A Mo - 0,84 et OC - 3 a d' incidence, le processus d'adptation converge 
en trois iterations vers une solution vraisemblable. 

an seul essai de decentrement de la maquette, sans doute trop important 
( _Z_ . . 57 Z), n'a pas permis de valider la methode dans ce genre de configu- 

ration. L. calcol -achfiaaclco. d. I’opfi refer 1“^" ^SSf f li. iion 

donne des solutions a caractSre oscillant, non utilisabJ.es pour 

correcte des deformations des parois. 
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Fig. 2 


ADAPTATION "3D" au centre veine 


Veine 2D 


parcis 


H+- B deformables 
laterales fixes // 


deplacement H , B = dep. sym. 4- dep. antisym. 




3D sur I'axe veine 


>y ♦ u = u,. +• 

x mesure interaction 
— *• parois 


U m. 

modeie 


Method e (exemple : sources =$> dep. svm .) 

*™sure_ UpW*^Cx-5)q(5)d5 = u H-t U B. 

-(=>«_ ♦ i nduft — u .( x )-y r K i ( x _ 5 ) q os)d5 

K p ,Kj fct (geometrie +• M 0 ) 
aveC /"KpCx-n) ft (n)dn = Kj(x) ft(t) , operateur lineaire 


Pour annuler I'interaction des parois u f C<) sur I'axe 

veine, on cree — Uj(x) par le de ptacement sv m. des 

parois Ail (x) = F[u D (x), ft(x), ] 

h p dx 2 

Demarche analogue pour dep. antisvm. avec 
tourbillons et u p (x) = u * 
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